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MASS-TO-LIGHT RATIOS OF ELLIPTICALS IN ACDM 

Napolitano N.R.^ and PN.S team^ 



Abstract. We use the mass-to-light gradients in early-type galaxies 
to infer the global dark matter fraction, fd =Md/Mt, for these sys- 
tems. We discuss implications about the total star formation effi- 
ciency in dark-matter halos and show that the trend of /d with mass 
produces virial mass-to-light ratios which are consistent with semi- 
analitical models. Preliminary kurtosis analysis of the quasi-constant 
M/L galaxies in Romanowsky et al. seems at odd with Dekel et al. 
simulations. 



^ 1 Introduction 

o 

. Recent finding of a sample of galaxies with a poor dark matter content within 
Qh' several effective radii (Romanowsky et al. 2003, R03 hereafter) seems to challenge 
Q the current paradigm of galaxy formation in ACDM cosmology where galaxies are 
?H ' seen to evolve in extended halos of dark matter with cuspy density profile (Navarro 
et al. 1997). Systems with a significant dark-matter component are expected to 
have mass-to-light ratios (M/L or T hereafter) deviating from the values typical 
of their stellar population. On the contrary R03 have shown that the kinematics 
of three intermediate luminosity galaxies can be modeled with a quasi-constant 
M/L consistent with the one expected from stellar population synthesis models. 
' This unexpected result have produced different interpretation either in the 
ACDM framework (Dekel et al. 2005 = DOS) or in MOND theory (Milgrom & 
Sanders 2003). In particular, D05 address very radial stellar orbits and projection 
effects in order to explain declining velocity dispersion profiles. 

On the other hand, Napolitano et al. (2005) (N05 hereafter) made predictions 
of gradients of mass-to-light ratios (V^T hereafter) in early-type galaxies and have 
ascertained that "quasi-constant M/L" are indeed expected within ACDM, albeit 
the R03 sample shows V^T which are too low and confiicting with acceptable star 
formation efficiency and baryon fraction. Low M/L gradients systems are a part 
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Fig. 1. Reduced kurtosis of NGC821 (full squares), NGC3379 (open squares), and 
NGC4494 (circles) compared to the average values from Dekel et al. 2005 (shaded boxes). 



of a generalised trend of V^T with stellar mass/luminosity (see Section 3). 
Here we show some preliminary results on the PNe's anisotropy of the R03 sample 
and compare them with the models from D05 and finally discuss whether the trend 
of the star formation eflaciency is consistent with virial M/Ls. 



2 Kurtosis profiles of PNe and anisotropy 

Simulations from D05 of galaxy merging of equal mass progenitors have shown 
that radial anisotropy naturally raises in the outskirts of merging remnants, where 
positive /i4 parameters should be measured. By converting /i4 to kurtosis using 
the approximated formula as in van der Marel & Franx (1993), it is possible to 
compare PN kurtosis to this predictions. In Fig. 1 we show the radial kurtosis 
profiles of galaxies NGC821, NGC3379 and NGC4494, obtained with the latest 
PN.S data reduction run compared with the quoted average values from Dekel et al. 
(represented as shaded boxes). It looks like the observed values, albeit consistent 
within the errorbars, are systematically lower than the prediction from DOS, which 
is already an indication of non strong radial anisotropy di per se. A direct estimate 
of orbital anisotropy needs data modeling: following Lokas & Mamon (2003) we 
have constrained the anisotropy parameter with Jeans analysis using 2nd and 4th 
velocity moments. Fig. 2 shows the; best-fit to the velocity dispersion and kurtosis 
for the case without dark-matter: best-fit parameters are T*_b = 12.8 ± 1.0 and 
/3 = —0.33 ± 0.12. Radial anisotropy {(3 > 0) is marginally allowed when forcing 
the models to a stellar mass-to-liglit ratio T*_b = 10 in order to ac;coinmo(late a 
dark halo. For this case halo concentration is c = 5, i.e. much lower then expected 
from N-body simulations (Bullock et al. 2001). 
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Fig. 2. Left: Velocity dispersion of NGC821 from long-slit (diamonds) and PNe (full 
squares) with best fit model without dark halo assuming isotropy (solid line) and (5 = 
—0.33 (dashed line). Right: reduced kurtosis of NGC821, symbols as center panel. 



3 Global dark-matter fraction and star formation efficiency 

In Fig. 3 (left) we compare the model V^T to the observed gradients of a sample 

of 21 early-type galaxies. As discussed in N05, V^T are considered as a two 
parameter family models (stellar, M*, and dark mass, or /d =Mi/M^ as 
adopted in Fig. 3). The dark-to- luminous fraction, /a, can be interpreted in 
terms of star formation efficiency, esp = -^*/-^bar, i-t^- the fraction of baryonic 
mass Mbar cooled in stars, assuming baryon conservation, such that esF = 4.9//d 
(see N05 for further details). 

In Fig. 3, the observed gradients do not seem to follow the simple increase of 
V^T with M« that is expected for a "universal" esF- This is made clearer in 
Fig. 3 (right) where the inferred esF (/d) are plotted against the galaxy masses: 
massive (brightest) sample is broadly distributed within the physical meaningful 
esF range marked by the dashed lines, while less massive galaxies (included the R03 
sample) have a very sharp distribution around a maximum esF(> 1) values which 
is unphysical. The mass scale marking this "dichotomy" is log(Mo/M0) ^11.2 
[Mq ~ 1.6 X 10^^Mq\. The overall trend seems to confirm other independent 
evidences that the star formation efficiency is a function of the halo masses with a 
minimum aroimd galaxies (Benson et al. 2000, Dckcl & Bornboi 2004, see also 
Marinoni & Hudson 2002, van den Bosch et al. 2003). The too high esF (too low 
/d) we have found around log(Mo/M0) 10.5 — 11.0 is still far to be understood 
and it remains a challenge for the theory. 

Driven by the results of Benson et al. (2000) and Dekel & Bornboi (2004), we can 
interpolate a smoothed /d function through the values plotted in Fig. 3 of the form 
/d = A{x - X(,f + B where x = \og{M^/MQY ■ Forcing the fit to a more physical 
minimum (we take fd = B = 4.9, i.e. esF = 1), the best fit (overplotted in Fig. 3, 
right panel) is found for xo = 10.9 and with A = 13.8 for a; < 11 and A = 68 for 



^The dependence of the esF on the mass can be obtained by the transformation egp = 4.9//d 
(see also N05). 
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Fig. 3. Left: model V^T (lines) and observed V^T (dots) from N05. Right: inferred 
/d = (esp) for individual galaxies (squaxes with errorbars) as in N05 with overplotted 
smoothed /a discussed in Section 3. See also text for details. 

X > 11. Is this result consistent with semi-analytical N-body simulations? 
Models from Benson et al. (2000) make explicit prediction of the M/L at the virial 
radius. Wc use the toy- model illustrated in N05 (in their Sect. 3) in order to 
compute the M/L at the virial radius produced by the /a function derived above. 
Indeed, by the definition of M/L: 

— = T* — = T* — =5.9 , (3.1) 

L M* eSFMbar esF 

where and Md are referred to the virial radius, T* is a constant stellar mass-to- 
light ratio (from population synthesis models) and we have used the cosmic baryon 
fraction fibar/^tot = 0.17. This shows that M/Ls at the virial radius depend on 
luminosity (and mass) through the quantities and esF- 

There arc indications that has only a weak dependence on (Trujillo et al. 
2004, we use their results in the following) which would leave esF as the main player 
of the virial M/L dependence on the mass. Virial M/Lb predicted according with 
Eq. 3.1 are compared with results from Benson et al. (2000) in Fig. 4, where 
we have rescaled their halo masses to stellar masses trough our /d. There is a 
substantial agreement with semi-analytical model prediction except in the low 
mass regime where our esF estimates arc based on poor data. But the overall 
behavior is reproduced. We also see that if we assume a stronger dependence 
of the T* on the mass (T» oc M^-^), the trend of our predicted virial M/L is 
inconsistent with Benson et al. (2000). 

4 Conclusions 

The main conclusions are: 1) for the moderate luminosity sample of R03, strong 
radial anisotropy does not seem to be an issue; 2) we have shown that the global 
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Fig. 4. Virial M/L derived from the /d smoothed function (solid line; long and 
short-dashed lines are uncertainties on the model) compared with the predictions from 
Benson et al. (2000) (data points with errorbars). Dotted blue line is the same as solid 
line but assuming T. cx M°-' instead of T* oc M° °^ from Trujillo et al. (2004). See also 
text for details. 

dark-matter fraction in ellipticals is a smoothed function of the stellar mass, con- 
sistently with the predictions of semi-analytical models; 3) the low-M/L gradients 
found by R03 and N05, for intermediate luminosity systems, imply too low dark- 
matter fraction challenging the standard ACDM, a problem that can be alleviated 
assuming low concentrated halos (R03, N05) or a lower as (Napolitano 2004). 
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